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^ ■ Abstract. 

We investigate the clustering properties of high-redshift galaxies within 
three competing scenarios for assigning luminous galaxies to dark matter 
O , halos from N-body simulations: a one galaxy per massive halo model, a 

^J ' quiescent star formation model, and a collisional starburst model. We 

^— N , compare these models to observations of Lyman-Break galaxies at 2; ~ 3. 

^ ' With current data and the simple statistic used here, one cannot rule 

out any of these models, but we see potential for finding distinguishing 
features using statistics that are sensitive to the tails of the distribution, 
and statistics based on the number of multiple galaxies per halo, which 

\j^ \ we explore in an ongoing study. 
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<3\ • !• Introduction 

Q^' Recent years have seen a new wealth of data on high-redshift galaxies. In par- 

' \ ticular, Steidel and collaborators have pioneered an efficient technique for iden- 

;_! ■ tifying a population of galaxies between redshifts of 2; ~ 2.5 — 3.5, and more 

'^ \ recently out to 2; ~ 4.5. These galaxies are known as the Lyman-break galaxies 

C^ ■ (LBGs) due to their color-selection at the Lyman-break. 

One of the first things noticed about these galaxies was that they are highly 
clustered, and several workers have investigated whether the clustering proper- 
ties of these galaxies could be used to distinguish between various cosmologies. 
C^ ' In Wechsler et al. (1998, hereafter W98) we compared the clustering of LBGs ob- 

served by Steidel et al. (1998, hereafter S98) to N-body simulations of 4 different 
CDM cosmologies. W98 and a number of other authors (e.g., Katz et al. 1998, 
Adelberger et al. 1998, hereafter A98) have found that the clustering properties 
of Lyman-break galaxies can be well matched by most reasonable cosmologies. 
Although it may be difficult to distinguish between cosmologies with this 
high-redshift data, there is still the hope that these galaxies may have much to 
teach us about galaxy formation. Most of the previous investigations of this 
subject have made simple assumptions about the relation between dark mat- 
ter halos and visible galaxies, and found that the data can be well-matched by 
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assuming a one-to-one correspondence between galaxies and the most massive 
halos. However, the detailed relationship is probably not this simple; in this pro- 
ceeding we investigate the effects of relaxing this assumption, and focus on the 
question of how the clustering properties of the LBGs can help us to understand 
their relation to dark matter halos and constrain models of galaxy formation. 
In this proceeding, we use N-body simulations to calculate the expected clus- 
tering properties of high-redshift galaxies under different assumptions for their 
halo occupation, motivated by several popular models of galaxy formation. Be- 
fore doing this, however, it is useful to investigate a simple model in order to 
gain understanding of how varying assumptions can affect the LEG clustering 
properties. 

There have been two general scenarios proposed for the nature of the LBGs. 
In one, known as 'Central Quiescent' (see, e.g., S98, A98, Governato et al. 1998, 
Baugh et al. 1998) the LBGs form only at the centers of the most massive galax- 
ies, which slowly accrete gas and form stars. In the other, known as 'Collisional 
Starburst' (e.g., Lowenthal et al. 1997, Somerville, Primack &: Faber 1999, here- 
after SPF, Kolatt et al. 1999) most of the LBGs are formed by small colliding 
galaxies that trigger a burst of star formation. Consider two simple models for 
the halo occupation in the above scenarios: first, for central quiescent LBGs, 

CO 

allow one LBG per halo above a mass threshold M^^, and second, for collisional 
starburst LBGs, let the average number of LBGs per halo of mass M > M^f^, 
scale like n{M) oc M^ , where S = 1.3 — this choice is motivated further in §4. 
In the first case, since LBGs are found only in the most massive halos, their 
clustering will be biased with respect to the underlying dark matter distribu- 
tion. In the second case, although collisions can be found in small-mass halos 
(typically M^f^ < M^^), they will preferentially be located within large halos, 
and the distribution of collisions should be also strongly biased. 

A quantitative estimate of this bias can be obtained using the model of Mo 
& White (1996), who give an expression for the halo bias bh = Ch/^DM- The 
appropriate bias factor for a population of galaxies within halos more massive 
than Mfnin is an average of bh{M,z) weighted by the abundance of halos as 
a function of mass, dN^/dM (e.g., as estimated by the Press-Schechter 1974 
approximation), and the average number of galaxies per halo n{M): 

1 C^ dNh 
bg{z, M > M„,„) = — — / -TTf{M, z)bh{M, z)n{M)dM, (1) 

where Ng{z) = J^ _ ■j^{M,z)n{M)dM. For the simple central quiescent sce- 
narios discussed above, n{M) = 1 and we obtain the standard expression for 
halo bias. We use expressions for dNh/dM and b^ supplied by Sheth &: Tormen 
(1999). 

In Figure la, we plot bg{Mmin) at z = 3 for both scenarios assuming the 
ACDM cosmology discussed below. Because high-mass halos are weighted more 
strongly in the collisional starburst model, we expect galaxies to be more biased 
for fixed Mmm- If we are free to set the mass threshold, we see that typical colli- 
sional starburst scenarios {M^f^ ~ 10^"'^^~^Mq ) and central quiescent scenarios 

i-^min ~ 10^^/i~^Mq ) can yield similar clustering properties. This is shown 
explicitly using N-body simulations by Kolatt et al. (1999). 
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Figure 1. Left: The bias parameter for halos and collisions as a 
function of Mmin- Right: The cumulative distribution of the error- 
weighted galaxy overdensity, for three LBG-halo models described in 
§4, compared with the data of A98. 



2. Data 

In the present analysis, we compare to the most recent published data of A98, 
who have found 376 galaxies in six 9' x 9' fields with redshifts between 2.5 and 
3.5. In analyzing this data, and 'observing' the simulations, we use a selection 
function fit to a histogram of all data, assume that 100% of galaxies with R > 
25.5 are photometric candidates at z ~ 3, and assume that ~ 45% of candidates 
have spectroscopic redshifts. 



Simulations 



The results we present here are based on N-body simulations of the GIF col- 
laboration, which were carried out at the Max-Planck-Institut fiir Astrophysik, 
Garching and the Edinburgh Parallel Computing Centre using codes kindly made 
available by the Virgo Supercomputing Consortium. The halo catalogues have 
been used for other purposes in, e.g., Kauffmann et al. 1999. Here we focus on 
one cosmology, an ACDM model with Q^n = 0.3 and h = 0.7. We use dark 
matter halos at z = 3 from a simulated box which is 141 /i~^Mpc on a side. We 
then test a series of models for populating these dark matter halos with visible 
galaxies, detailed in the following section. 

Once the halos have been populated with galaxies, we observe these galaxies 
exactly as the data is observed, making the same assumptions about the selection 
function mentioned above. We normalize each model to get the correct number 
density of observed galaxies. This normalization is different for each type of 
model, and is discussed further below. 



4. Models for populating halos 

(1) One LBG per massive halo: This simple model for populating halos 
with Lyman-Break galaxies is motivated by the assumption that the luminosity 
of these galaxies is directly proportional to mass, and has been considered by 
several authors (e.g., W98, Jing & Suto 1998). In this model the number density 
is fixed to match the data by choosing an appropriate Mmm- 

(2) Semi-analytic models with quiescent star formation: Since 
clearly the model discussed above is somewhat simplistic, we consider a more 
accurate way to fill halos with quiescently star- forming galaxies, using the semi- 
analytic models of SPF, which give the probability of observing galaxies of var- 
ious luminosities given the mass of the halo. Here we only consider one model, 
referred to as SFR-M in SPF; see Wechsler et al. (1999, hereafter W99) for 
consideration of a wider range of models. In this class of models, the number 
density of objects is normalized using an adjustable dust parameter described in 
SPF &; W99. The dust parameter we use is slightly greater than that estimated 
by Steidel et al. (1999) for these galaxies. 

(3) A toy model for collision-driven starbursts: In the collisional 
starburst picture, most of the visible high-redshift galaxies are starbursting due 
to recent collisions. We have analyzed a very high-resolution (ART) N-body 
simulation of the same ACDM model within a 30 /i~^Mpc box (Kravstov, Klypin 
& Khokhlov 1997) where halo substructure and mergers have been identified (for 
details see Bullock et al. 1999, Kolatt et al. 1999). Using collisions identified in 
this simulation, we have found that the number of collisions per host halo, ricou, 
follows a power law with host halo mass Mhostj 't-coU cx^ -^i^st' with S ~ 1.3. 
For further details on how this model is fit from simulations and how it might 
be explained, see W99. We use this power-law model to populate the GIF dark 
matter halos with collisional starburst galaxies. We take the minimum mass for 
a halo to host a collision to be Mmin = lO^^/i~^M0; the number density for this 
model is set by adjusting the constant in front of the power law. 

5. Comparison of models with data 

Here we just consider one statistic for comparing our models with the current 
data; for a more thorough analysis see W99. A standard statistic for measuring 
the clustering of a population is the overdensity in some region; here we consider 
the overdensity in cells that are 5z = 0.04 in redshift and 9' x 9' in angle. In 
order to minimize the contribution of cells where the selection function is very 
small, we instead look an error-weighted overdensity per pixel, defined as: 

D, = 5N,/ai = 5N,/Nj^, (2) 

where Ni is the number per pixel and 5Ni = {Ni — N)/N . 

For each model, we compare the cumulative distribution of this statistic 
with that of the data (Figure lb), and there is some indication that the colli- 
sional model has an excess of both overdense and underdense cells relative to 
the data. A detailed Kolmogorov-Smirnov test for each model, which gives the 
probabilities that the data and the model came from the same underlying dis- 
tribution, yields: 



model (1): Pks = 0.15 ± 0.06 

model (2): Pks = 0.23 ± 0.10 

model (3): Prs = 0.10 ± 0.04 
In this analysis, there are 144 data pixels and 720 simulation pixels; we do the 
assignment of galaxies to halos and 'observation' of LBGs 10 times for each 
model to get an error estimate. 

6. Discussion 

The results above show that none of the three models can be rejected with this 
statistic and the current data, but there are indications that the different models 
may provide distinguishable clustering properties. In W99, we consider both a 
wider range of models and more discriminating statistics, including the number 
of close pairs and a Kuiper test, sensitive to distribution tails. 
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